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We continue the systematic computation of Landau gauge gluon and ghost propagators of SU(2) 
gluodynamics using a sequence of increasing lattice sizes L 4 up to L = 112 with corresponding 
j3 -values chosen to keep the linear physical size fl(j3)L ~ 9.6 fm fixed. To extremize the Lan- 
dau gauge functional we employ simulated annealing combined with subsequent overrelaxation. 
Renormalizing the propagators at momentum ji = 2.2 GeV we observe quite strong lattice arti- 
facts for the gluon propagator as well as for the ghost dressing function within the momentum 
region q < 1 .0 GeV. The dependence on the lattice spacing for the gluon propagator at lowest 
accessible physical momentum values does not yet allow a simple extrapolation to the continuum 
limit. On the contrary, the running coupling derived from the bare dressing functions seems less 
affected by lattice artifacts. 
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Table 1: Parameters of the main simulations considered here. See also the text. 




Figure 1: Left: the unrenormalized gluon propagator D(q 2 ) and its fits according to Eq. (1); the data points 
drawn at q 2 = 0.001 represent the zero-momentum gluon propagator D(Q) values. Right: the ghost dressing 
function J(q 2 )- Both are shown for approximately equal physical volume ( (aL) 4 ~ (9.6fm) 4 ) but different 
j3 -values, i.e. discretization scales a. 



The aim of the reported study is to continue the systematic investigation of SU (2) gluon and 
ghost Landau gauge propagators on large lattices [1] in order to receive information from first 
principles on the behavior of these propagators and of the running coupling in the continuum limit 
for all momenta q including the infrared (IR) region. As in previous investigations [2], for solving 
the Gribov problem we assume the Landau gauge functional to be driven as close as possible 
to its global extremum. Employing the standard Wilson plaquette action we have studied gluon 
(D) and ghost (G) propagators for lattice sizes L 4 with run parameters collected in Table |TJ The 
gluon propagator was computed using Nmc(P,L) independent Monte Carlo (MC) configurations 
generated with the given set of parameters, while the ghost propagator was calculated only on a 
subset of N g h ost (p,L) MC configurations. We used the same procedures for Landau gauge fixing 
and computation of propagators as described in [1], namely, we employed very long simulated 
annealing (SA) runs followed by overrelaxation (OR) to obtain gauge copies with a gauge fixing 
functional close to its global extremum for each MC configuration. In Fig. [I] we show the results 
for the bare gluon propagator (left) and the bare ghost dressing function (right) for fixed physical 
volume (aL) 4 but varying lattice scale a. For the gluon propagator we have drawn also curves 
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Table 2: Results of the 6-parameter fits of the unrenormalized gluon propagator for various /3 corresponding 
to lattices sizes 56 4 , 80 4 , 96 4 and 1 12 4 . 



obtained from fits with the 6-parameter formula proposed in [3] 

q*+A 2 q 2 +B 

We found the resulting fit curves nicely to capture the IR turnover of the gluon propagator. The 
% 2 /dof values are close to unity in most cases (see Table ^). To obtain the renormalized gluon 
propagator D ren (q,n) = Z(a,fi)D(q,a)) we apply the normalization condition D ren (/x,/x) = 1/jU . 
Since the fit formula Eq. (|l|) nicely works throughout the whole momentum region we can use it to 
carry out the renormalization at any jU. For ju = 2.2 GeV the renormalization factors Z are collected 
in Table []]. Their values do not vary strongly which means that the bare dressing functions for the 
different /3-values approximately overlap at the given jU-value. From Fig. [I] we can then conclude 
also for the renormalized gluon propagator D re „(q ) and the renormalized ghost dressing function 
J r en{q 2 ) found for various lattice spacings aQ3) to be compatible with the so-called decoupling 
solution of Dyson-Schwinger or functional renormalization group equations (see [4]). The numer- 
ical values, however, of D ren (q 2 ) and J ren {q 2 ) in the limit q — > appear to be j8- or a-dependent. 
From such plots one can see that the convergence of the renormalized lattice propagators / dressing 
functions in the deep IR momentum range to the respective continuum counterpart, that should be 
observed for decreasing is rather slow. For their direct numerical study near the continuum 
limit one has to use rather large j8 -values which consequently requires simulations on unrealisti- 
cally huge lattices, which are not accessible today even on most powerful parallel supercomputers. 
Instead, we can try to make contact with the continuum limit by extrapolating D ren (q 2 ,a) to the 
zero-a limit as done e.g. in [5] for SU (3) and non-zero temperature. In Fig. || (left) we plot the 
a-dependence of lattice D rm (q 2 ) for several selected values of q 2 . We see the lower the momentum 
is the less well-defind the convergence for a — > becomes. For getting reliable numerical values 
of SU (2) gluon and ghost propagators in the continuum limit more work is needed. Although 
the ghost dressing function J(q 2 ) has been computed only for a subset of MC configurations (see 
Table ||), it provided useful quantitative information, see Fig. [I] (right). Even a single MC con- 
figuration as for L = 96 already seems to yield a first estimate (a fast decrease of the statistical 
fluctuations of J(q 2 ) with increasing L was first observed for the SU (3) case in [6]). Our analysis 
shows in the deep IR region that, while D ren {q 2 ) increases with j8, J re n{q 2 ) decreases. More details 
will be published elsewhere. 
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Figure 2: Left: Renormalized gluon propagator D ren (q 2 ) at jx = 2.2 GeV versus lattice spacing a for various 
physical momenta q 2 (values indicated in units GeV 2 ). Right: Running coupling for /3 = 2.3, /3 = 2.4 and 
J3 = 2.45 




Figure 3: Check of systematic errors. Left: The unrenormalized D(q 2 ) computed for j3 = 2.3 and various 
L. Right: D(q 2 ) for 2 different SA schedules at J3 = 2.4 and L = 80. 



We have checked whether the differences of propagator values in the deep IR could be com- 
pensated by other systematic effects. From Fig. [| (left) one can see that finite- volume effects are 
small if the linear physical size is a(j3)L ~ 9.6 fm or even larger. What concerns Gribov copy 
artifacts at /3 = 2.4 and L = 80 we have compared the results of two sets of SA+OR gauge fixing 
simulations: (i) one gauge copy fixing with 9600 SA sweeps ("SA1 schedule") (Nmc = 314) and 
(ii) "best of two copies" gauge fixing with 12000 SA sweeps each ("SA2 schedule") (Nmc = 187). 
For details see, e.g., Ref. [7]. SA1 and SA2 results obtained for unrenormalized gluon propagators 
are plotted in Fig. ||| (right). We have found that the differences between these cases are much 
smaller than the magnitude of Gribov copy effects measured in [1] as difference between results 
of one-copy SA+OR and one-copy OR gauge-fixing procedures. Our analysis shows that further 
"improvement" of SA schedules could not change D ren (q 2 ) essentially and hence noticeable dif- 
ferences of D ren (q 2 ) values in the deep IR region found for different j3 -values certainly cannot be 
accounted for by the Gribov copy effect. 
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With the bare gluon (Z(q 2 )) and ghost (J(q 2 )) dressing functions at hand one can easily com- 
pute the running coupling in the minimal MOM scheme [8], 

OC s (q 2 ) = ^J 2 (q 2 )Z(q 2 ). 

The dependence of the resulting curves on j8 or a turns out to be rather weak even in the deep 
IR momentum region (see Fig. || (right)), i.e. the lattice artifacts of the gluon and ghost dressing 
functions cancel each other to some extent. 

We conclude that naive multiplicative renormalizability for the SU (2) Landau gauge gluon and 
ghost propagators gets violated in the deep IR region. Due to the slow convergence of gluon and 
ghost renormalized propagators their continuum counterparts may strongly differ in the deep IR 
momentum region from what we have obtained here in lattice simulations with admissible values 
of /3 = 4/#q. At the same time, the physically important renorm-invariant minimal MOM-scheme 
running coupling a s (q 2 ) seems to reach a continuum behavior much earlier. 

IB thanks Prof. A. A. Slavnov for a useful discussion of the results. Simulations have been 
done on the MVS100K supercomputer of the Joint Supercomputer Centre (JSCC, Moscow). 
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